Exotic phenomena in doped quantum magnets 

o. 

^ : D. Poilblanc\ M. Mambrini^ A. Lauchli^ F. Mila^ 

^ Laboratoire de Physique Theorique, CNRS & Universite de Toulouse, 118 rte de 
Narbonne, 31062 Toulouse, France 

^ Institut Romand de Recherche Numerique en Physique des Materiaux (IRRMA), 
PPH-Ecublens, CH-1015 Lausanne, Switzerland 
QJ ■ ^ Institute of Theoretical Physics, Ecole Polytechnique Federale de Lausanne, BSP 

^ \ 720, CH-1015 Lausanne, Switzerland 

^ . E-mail: didier.poilblanc@irsamc.ups-tlse.fr 

^ I Abstract. We investigate the properties of the two-dimensional frustrated quantum 

^ ■ antiferromagnet on the square lattice, especially at infinitesimal doping. We find that 

^ . next nearest neighbor (N.N.) J2 and next-next N.N. J3 interactions together destroy 

O ■ the antiferromagnetic long range order and stabilize a quantum disordered valence 

bond crystalline plaquette phase. A static vacancy or a dynamic hole doped into 
this phase liberates a spinon. From the profile of the spinon wavefunction around the 

■ (static) vacancy we identify an intermediate behavior between complete deconfinement 
. (behavior seen in the kagome lattice) and strong confinement (behavior seen in the 

QO ' checkerboard lattice) with the emergence of two length scales, a spinon confinement 

■ length larger than the magnetic correlation length. When a finite hole hopping is 
introduced, this behavior translates into an extended (mobile) spinon-holon boundstate 

\^ ■ with a very small quasiparticle weight. These features provide clear evidence for a 

■ nearby "deconfined critical point" in a doped microscopic model. Finally, we give 
' arguments in favor of superconducting properties of the doped plaquette phase. 
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Figure 1. (Color on-line) Exemples of simple VBC states, (a) columnar phase and 
(b) plaquette phase, where 2 (4) spins are paired up in dimer- (plaquette-) singlets. In 
both cases the ground state is 4-fold degenerate. The J2 and J3 couplings (depicted 
on figure (b)) are shown to stabilize the plaquette phase. 



1. Introduction: the plaquette valence bond crystal 

Magnetic frustration is believed to be the major tool to drive a two-dimensional (2D) 
quantum antiferromagnet (AF) into exotic quantum disordered SU(2)-symmetric phases 
such as the spin liquid (SL) state characterized by the absence of ordering of any kind, 
and possibly observed in the 2D kagome lattice [Ij. The Valence Bond Crystal (VBC) 
which, in contrast, breaks lattice symmetry (see Fig.n]for pictures of such states), seems 
to be a serious alternative in some other frustrated quantum magnets as suggested 
by robust field theoretical arguments early numerical computations of frustrated 
quantum AF on the square lattice with diagonal J2 bonds and in the 2D checkerboard 
lattice [4] (with diagonal bonds only on half of the plaquettes). 

The "deconfined critical point" (DCP) was recently proposed to describe a new 
class of quantum criticality characterizing the AF to VBC transition Here we argue 
that the two-dimensional (2D) spin-1/2 AF J1-J2-J3 Heisenberg model (Hj on the square 
lattice defined by 

H = ^JijSi'Sj (1) 

{ij) 

where the Jij exchange parameters are limited to first (Ji), second (J2) and third 
(J3) N.N. AF couplings is a strong candidate for exhibiting such a transition. Recent 
investigations of this model along the pure J3 axis [3 [Hj or for both J2 and J3 finite [9j 
provide an accumulation of evidence in favor of the existence of a VBC phase with 
strong plaquette correlations similar to the cartoon of Fig. Q^b). 
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Figure 2. (Color on-line) Classical phase diagram for the J1-J2-J3 model. Second 
order (discontinuous) transitions are indicated by dashed (solid) lines (see Ref. 6 for 
more details). The shaded (blue online) region shows the approximate location of the 
minimum of the impurity spectral weight Z (0.79 < Z < 0.84) in the quantum version. 

Extensive exact diagonalisations of the model within its full Hilbert space and 
within a restricted Hilbert space of N.N. (SU(2)) dimer configurations reported in Ref. [9j 
strongly support the existence of a plaquette VBC phase along the J2 + J3 = 1/2 hne 
of parameter space (see Fig. Ej). We shall provide here a selection of these data to 
illustrate this point. Dimer-dimer correlations (not shown) reveal a strong signal of 
VBC order of some kind. Summing up these spatial correlations with appropriate phase 
factors provides quantitative estimates of both the "generic" VBC and the "specific" 
columnar VBC structure factors (normalized to give the squares of the corresponding 
order parameters). Such estimates are displayed in Fig.O the upper line is obtained by 
an extrapolation to infinite size of the dimer structure factor constructed by summing up 
the correlations of the horizontal dimers (only) with an alternating sign along columns 
(in order to provide a signal for all types of VBC phase). Clearly, a large VBC order 
parameter is seen with a value close to the one of a pure plaquette phase (dotted line), 
at leat close to the pure J3 axis. The data sets at the bottom are obtained by summing 
up the dimer correlations with opposite sign for vertical and horizontal dimers (also 
omitting the short distance contributions) to filter the signal of columnar order only. 
The small values of the latter data strongly suggest that the VBC order is a plaquette 
order rather than a columnar order. Note, however, that the VBC signal weakens 
significantly when approaching the pure J2 axis. 




Figure 3. (Color on-line). VBC (extrapolation) and columnar (for 32 and 40 sites 
clusters) structure factors along a simple cut J2 + J3 = 1/2 in parameter space (see 
text). The dashed lines correspond to the value of both structure factors expected in 
pure columnar dimer and plaquette states. From Ref.j2|. 



2. Doping the plaquette VBC 

We shall now test the properties of the plaquette VBC by removing an electron at a given 
site (e.g. by chemical substitution with an inert atom) or, as in ARPES experiments, 
in a Bloch state of given momentum. This process naturally liberates a spinon, i.e. a 
S=l/2 polarisation in the vicinity of the empty site (holon). A new length scale ^coni 
corresponding to the average distance between vacancy and spinon emerges naturally in 
a VBC phase and is to be identified with the correlation length over which VBC order 
sets in. Interestingly, it has been predicted that, in the vicinity of the DCP, confinement 
occurs on a much larger length scale ^conf than the spin-spin correlation length ^af [5j . 

In Fig. HI we compare the decay of the spin-spin correlation with distance to the one 
of the spinon " cloud" away from the vacancy which both can be fitted usig exponential 
forms, hence enabling to extract the corresponding length scales. The obtained very 
short magnetic correlation length (^af, below one lattice spacing, is to be contrasted with 
the strikingly large confinement length ^conf typically ranging from 2 to 6 lattice spacings. 
This is to be compared with two other behaviors observed in the kagome and in the 
checkerboard lattices [lOj respectively: in the checkerboard lattice (whichalso exhibits 
a plaquette phase but a larger spin gap ^) very short-ranged spin-spin correlations are 
seen while the spinon remains almost entirely confined on the N.N. sites of the vacancy. 
In contrast, on the kagome lattice, the spin-1/2 delocalizes on the whole lattice, a clear 
signature of deconfinement which is consistent with a SL ground state. 

The quasiparticle weight Z can be deduced from the computation of the single hole 
Green function. Fig. [21 shows that it is clearly reduced in the vicinity of the J2+ J3 = 
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Figure 4. (Color on-line) Modulus of the spin-spin correlation (black) and of the spin 
polarization in the vicinity of the vacancy (red) (summed up on equivalent sites) versus 
distance for both J2/J1 = 0.3 and J3/J1 = 0.2 (a) and J2/J1 = 0.1 and J3/J1 = 0.4 (b) 
corresponding to points A and B in the phase diagram of Fig. [21 Fits using exponential 
forms are shown in dashed lines. The areas of the dots are proportional to the number 
of equivalent sites (entering in the fits). Dark and light symbols correspond to positive 
and negative values respectively. From ref. pUj . 



Jl/2 line where evidence for a plaquette VBC phase was found. This behavior is in 
agreement with the above real-space picture which naturally implies that the extended 
spinon wavefunction of size ^conf has a reduced overlap Z with the bare wavefunction (of 
extension <^af)- When the vacancy is given some kinetic energy (through a N.N. hopping 
as in the well-known t-J model), the Z factor of the hole drops further to very small 
values [ini showing that the hole motion strongly suppresses the effect of the remaining 
spinon-holon VBC string potential [llj. Typical spectral functions along the pure J3 
line are shown in Fig. El Note that a complete deconfinement was shown for a mobile 
hole on the kagome lattice fT^ as evidenced by a fully incoherent spectral weight at low 
energies. 

Lastly, we would liketo comment on the possibility of superconducting pairing 
upon doping VBC's. In fact, pairing driven by a kinetic energy gain of the Cooper 
pair was discovered in the robust plaquette phase of the checkerboard lattice [13j. On 
general grounds, one would also expect hole pairing in the doped J1-J2-J3 quantum 
antiferromagnet due to the long distance confining string of the VBC. However, short 
range effects (as in the checkerboard lattice) could boost the weak pairing interaction. 
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Figure 5. (Color on-line). Single hole spectral function (at momentum (7r,0)) for 
various values of the magnetic frustration and a realistic hole hopping of t/Ji = 2.5. 
Note that for J3/J1 = 0.5 the main low energy peak contains several poles in contrast 
to the J3 = 0.05 case where a single pole is present. 



3. Conclusions 

To conclude, the confinement of a spinon liberated by introducing a vacant site or 
a mobile hole has been studied in the J1-J2-J3 model where a plaquette VBC phase 
has been identified in some extended region of the phase diagram. In this region of 
large frustration, we have identified a new length scale related to the confinement of 
the spinon. Its large value compared to the spin-spin correlation length supports the 
field-theoretic "deconfined critical point scenario" [3] for the Neel-VBC transition. 
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